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We extract a set of values for the Gross-Llewellyn Smith 
sum rule at different values of 4-momentum transfer squared 
(Q 2 ), by combining revised CCFR neutrino data with data 
from other neutrino deep-inelastic scattering experiments for 
1 < Q 2 < 15 GeV 2 /c 2 . A comparison with the order a 3 theo¬ 
retical predictions yields a determination of a s at the scale of 
the Z-boson mass of 0.114±;oi2 . This measurement provides 
a new and useful test of perturbative QCD at low Q 2 , because 
of the low uncertainties in the higher order calculations. 

PAC Numbers: 12.38.Qk 11.55.Hx 13.15.+g 24.85.+P 


The Gross-Llewellyn Smith (GLS) sum rule [Q predicts 
the integral Jg(xF 3 ) — , where xF 3 {x,Q 2 ) is the non¬ 
singlet structure function measured in neutrino-nucleon 
(yN) scattering. In the naive quark parton model, the 
value of this integral should be three, the number of va¬ 
lence quarks in the nucleon. In perturbative Quantum 
Chromodynamics (pQCD), this integral is a function of 
a s (Q 2 ), the strong coupling constant. 

The GLS integral is one of the few physical quanti¬ 
ties which has been calculated to next-to-next-to-leading- 
order (NNLO) of perturbative QCD Q, and there are 
estimates of the next order term j3] (i.e. 0(a 4 )). In ad¬ 
dition, there is a non-perturbative higher-twist contribu¬ 
tion, proportional to 1 /Q 2 . This yields the GLS integral 
as a function of a s , of the form: 


GLS = 3 
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where a(rif) and b(rif ) are functions Q of the number 
of quark flavors accessible at a given Q 2 . The higher- 
twist correction term A HT is predicted to be significant 
in some models Q, while others predict a negligibly 
small correction term. We take A HT as half the largest 
model prediction, with errors which cover the full range 


(ART = 0.15±.15 GeV 2 ). 

The size and Q 2 -variation of the GLS integral is a ro¬ 
bust prediction in pQCD. The NNLO calculation has 
been shown to be largely independent of renormaliza¬ 
tion scheme j|] , and xF 3 is inherently independent of the 
gluon distribution. The number of orders to which the 
integral has been calculated ensures an accurate pertur¬ 
bative calculation in spite of the large value of a s at low 

Q 2 - 

An earlier measurement of the GLS integral has been 
published by the CCFR collaboration |)j. That analysis 
used a leading-order(LO) QCD-based fit to extrapolate 
all data to Q 2 = 3 GeV 2 , fitted the extrapolated data to a 
single function over all x, and numerically integrated that 
function. This was confirmed by a LO global fit analysis 
of the same data ©• However, these approaches cannot 
make full use of the accuracy of the NNLO calculation 
shown above, since they depend on LO pQCD for extrap¬ 
olation. Also, the previous CCFR analysis did not correct 
for quark mass thresholds ||, target mass |ll| or higher- 
twist 00 effects. These corrections are important at 
the effective mean Q 2 of the result (Q 2 ~ 3 GeV 2 ). 

This paper describes a new GLS analysis, which uses 
revised CCFR xF 3 data together with data from ear¬ 
lier neutrino-scattering experiments. By combining data 
sets, we expand the kinematic region to measure / xF 3 ^ 
for 1 < Q 2 < 15 GeV 2 without any extrapolation in Q 2 . 
This technique thus allows us to consistently use the fun¬ 
damental NNLO prediction shown in equation (1). 

The CCFR data were collected at Fermilab in exper¬ 
iments E744 and E770, which ran in 1985 and 1987-8 
respectively. The experiments observed neutrino scatter¬ 
ing in an iron calorimeter |l2[ . The calorimeter and muon 
spectrometer were calibrated using a test beam |li|. New 
structure functions (SFs) from this data [Q were pub¬ 
lished in 1997. These SFs had a number of improvements 
compared to the SFs used in the previous GLS measure¬ 
ment |9]. Improvements include a revised energy cali- 
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bration based directly on test-beam data, an improved 
calculation of radiative corrections |h| , and the removal 
of and correction for two-muon events (vN —> fi~X) 
from the data sample. Previously, the two-muon events 
introduced a small ambiguity at low x between neutrino- 
induced and anti-neutrino-induced events, which is par¬ 
ticularly important to the GLS integral. 

This analysis further improves the CCFR structure 
functions [ p~| at low x by improving the acceptance and 
smearing corrections. These corrections, which require 
a cross-section model, now incorporate measurements of 
the strange sea 116 and a more accurate parameteriza¬ 
tion of the parton distributions. These procedures create 
a new SF set jT^] with reduced uncertainty at low x. 

We also expand our kinematic region at high x by 
using xF 3 data from other i/-N experiments, namely: 


WA25 @, WA59 pi SKAT FNAL-E180 and 
BEBC-Gargamelle |22(|. These were each normalized to 
CCFR in the regions of overlap, and the WA25 data were 
corrected at high x for nuclear differences J|iJ in the tar¬ 
gets. 

The GLS integral is evaluated numerically using the 
combined xF 3 data in bins of x and Q 2 . The integral 
over x is evaluated separately for each Q 2 bin. At very 
low x we must extrapolate below the CCFR kinematic 
limit, while at high x we use other experiments’ data, 
and interpolate as necessary within the large bins. In 
each case, we vary the forms of the interpolations and 
extrapolations, and use the differences in the integral as 
estimates of the systematic uncertainties in the proce¬ 
dures. 

The CCFR data have a minimum x of roughly (x = 
0.002 x Q 2 ). To extrapolate below this, we fit a power law 
( Ax B ) to all points with x < 0.1. The power-law form is 
suggested by Regge theory [[m|, which predicts a shape 
of a; 0,5 . To test this assumption, we made an alternate fit 
of the form Ccc 0 ' 5 , using the difference as an independent 
systematic uncertainty. This systematic error becomes 
large at Q 2 above 5 GeV 2 . 

For x > 0.5, there are also few data points. Most of 
the data here come from BEBC and SKAT which quote 
only two points for the range 0.5 < x < 1.0. Here xF 3 is 
steeply falling and thus the precise shape is important for 
integration. Again, to estimate the contribution and er¬ 
ror we use various assumptions. For the central value, we 
use the principle that at high x, the shape of xF 3 should 
be the same as F 2 , since the sea quarks are negligible at 
high x. Electron scattering experiments at SLAC have 
precisely measured F- 2 in this region j2f|. These data are 
corrected for nuclear effects ^3| and differences between 
eN and vN scattering Q. The corrected F 2 data by 
itself give the same result as interpolating the xF 3 data 
with the F 2 shape. 

However, the SLAC data have small resonance peaks 
which may be different from neutrino resonances. To 
estimate the systematic error, we take the difference be¬ 


tween two power-law fits to the xF 3 data, using the forms 
D( 1 — x) E and F( 1 — a;) 3 . These bracket the SLAC fit 
and serve as the limits of reasonable interpolation forms. 
Note that resonance behavior at low Q 2 coupled with 
approximate scaling Q leads to a predicted form of 
(1-x) 3 . 

Fig. [t] shows the combined xF 3 data on a log x scale 
in four low-Q 2 regions, along with a line representing the 
power law fit ( Ax B ) for x < 0.1 and the y 2 for the fits. 
Fig. [ 2 ] shows the same data on a linear x scale to highlight 
the high-x region. 

Following the procedure of the BEBC collaboration 
p2f , we add the quasi-elastic contribution and correct the 
GLS integral for target mass effects. This is necessary to 
be consistent with theoretical prediction of higher-twist 
contributions || to the GLS integral. Table | shows the 
exact Q 2 ranges of each bin and the contributions to 
f xF 3 — for the different regions of x. 

The systematic uncertainties are divided into three 
classes. The first includes calibration, normalization, and 
other purely experimental issues. The second is uncer¬ 
tainty in the integration of experimental xF 3 , estimated 
by varying the assumed functional forms as described 
above. The third class includes uncertainties in the theo¬ 
retical prediction of the GLS integral itself. A summary 
of all the uncertainties is shown in Table ||. 

The dominant experimental systematic uncertainties 
are in the normalization of xF 3 , which comes from the 
total neutrino and anti-neutrino cross-sections (oy and 
ay). The absolute oy is not measured by CCFR, so 
we use the world average |27],|8) (a V /E v = 0.677 ± 
.007 x 10 _38 cm 2 /GeV). The ratio oy/cr v is measured 
by CCFR, and combined with the world average [E8 
yields ay/ay = 0.499 ± .007. Other experimental uncer¬ 
tainties include the energy scale calibration of the detec¬ 
tor and the effects of charm production on the measured 
structure functions. 

Additionally, there is a small uncertainty in the re¬ 
vised calculation of acceptance and smearing corrections. 
These corrections depend on a parameterization of the 
SFs. Variations of the functional form of the parameter¬ 
ization were used to estimate the systematic error. 

The dominant theoretical uncertainty is the error 
on the higher-twist correction (A HT/Q 2 ). Braun and 
Kolesnichenko Q use three models which predict a cor¬ 
rection term A FIT between 0.16 and 0.29 GeV 2 . Other 
models, such as bag models pj and a recent NNLO analy¬ 
sis ]g|] using a renormalon ]?]] approach, predict a negligi¬ 
ble correction term (A HT < 0.02 GeV 2 ). For our central 
value, we take AHT = 0.15 ± .15 GeV 2 , thus covering all 
three predictions. The nuclear effects of the target are 
predicted to be small |jo| (—0.01 /Q 2 for iron). We also 
use estimates of the renormalization scheme dependence 
0 and the order af term in the pQCD expansion || as 
uncertainties in the perturbative calculation. 
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To extract a single value for A-—!-, we combine the 
° MS 

measured values of the GLS integral in each Q 2 with 
the uncorrelated systematic errors, including the accep¬ 
tance model error and the high-a; and low-:r fitting errors. 
These points are fitted to the NNLO pQCD function 
and higher-twist term, shown in Eq. (.l[). The prediction 
includes quark mass thresholds using the procedure of 
Chyla and Kataev ||. The other systematic error sources 
are fully correlated in Q 2 , and are applied by shifting all 
GLS values by the uncertainty from that source and re¬ 
doing the fit. The difference between the shifted and 
unshifted fit result represents the uncertainty in a s from 
that systematic error. 

The best fit to the measured GLS integral as a function 
of Q 2 is for a value of A^L = 165 MeV. Evolving to M§ 
and 3 GeV 2 at NNLO, this corresponds to: 

a s (M|) = 0.114 ±;Qog (. stat ) ±;°°9 (syst) ± .005 ithry) (2) 
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a s (3 GeV 2 ) = 0.28 ± .035(sfof) ± .05 {syst) ± ( thry ) (3) [11] 

[ 12 ] 


If the higher twist models of [5-7] are used, the central 

values become a s (M§) = 0.118 and a s (3 GeV 2 ) = 0.31. 

Table III shows the results of fitting for ALjL at each 

° MS 

Q2 value as a consistency check. In all cases the small 
target mass and quasi-elastic corrections are included, 
and roughly cancel. 

In conclusion, the GLS sum rule allows a precise mea¬ 
surement of a s at low Q 2 . An independent measurement 
of a s from the CCFR calculation [Q used the slope of 
global NLO fits to xF 3 and F 2 for 15 < Q 2 < 125 GeV 2 , 
and found a s (M = 0.119 ± .004. Three orders of mag¬ 
nitude higher in scale ( Q 2 = M§ = 8315 GeV 2 ), elec- 
troweak fits to LEP data j3lj found a s (M§) = 0.124 ± 
.004 ± .002 based on the parameter Re. 

The GLS result is consistent with other measurements, 
showing the power of pQCD across a very wide range of 
scales. An inconsistency in results might have indicated a 
need for higher order calculations in other measurements 
(i.e . NNLO) or the presence of other theoretical effects 
[p2| which scale with Q 2 . The GLS result can be im¬ 
proved in the future by additional data from the NuTeV 
experiment at Fermilab and by better understanding of 
higher-twist effects. 
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FIG. 1. XF 3 as a function of x at the four lowest Q 2 val¬ 
ues, with x on a log scale. The area under the points thus 
represents J xFs^. The curve is a power law (Ax B ) fit to the 
x < 0.1 points, which is used to calculate the integral in the 
shaded region. 



FIG. 2. XF 3 as a function of x at the four lowest Q 2 values, 
with a: on a linear scale to show the high-x data. The shaded 
region shows the fit using the shape from SLAC F 2 data. 
The other lines are the power law fits (D( 1 — x) E above and 
F( 1 — x) 3 below) used to estimate systematic error. 


TABLE I. The contributions to the GLS integral from dif¬ 
ferent regions of x and the quasi-elastic peak added at x = 1, 
shown as a function of Q 2 (in GeV 2 ). For high and low x, 
it shows the estimated uncertainties due to the model choice 
only. The values include the target mass corrections. 


Q 2 

Jo 02 Fzdx 

fA^dx 

/.s' F3dx 

qElas 

1 . 0 - 1.6 

0.376 ± .082 

1.730 

0.183 ± .073 

0.103 

1 . 6 - 2.5 

0.523 ± .002 

1.843 

0.091 ± .026 

0.033 

2 . 5 - 4.0 

0.558 ± .026 

1.889 

0.092 ± .020 

0.009 

4.0- 6.3 

0.700 ±.137 

1.991 

0.084 ± .016 

0.002 

6.3-10.0 

0.748 ± .139 

2.004 

0.064 ± .008 

0.0004 

10.0-15.5 

0.718 ± .113 

2.007 

0.073 ± .003 

0.0001 


TABLE II. Uncertainties in a s (Mf). Errors which are un¬ 
correlated in Q 2 are marked with a *. Other sources are fully 
correlated in Q 2 . 


Source 

Error 

Statistical 

( + .005 \ 

v —.006 / 

cr„ Normalization 

( +.003\ 

V —.005 / 

a v 1 av Ratio 

( + .005 \ 

V —.006 / 

Energy Calibration 

( +.002\ 

^ —.003 / 

Charm Production 

±.0005 

Acceptance model* 

±.002 

Total Experimental Error 

( +.006\ 

^ —.008 / 

High-x fitting* 

±.003 

Low-x fitting* 

±.002 

Total Model Error 

±.004 

Combined Systematic Error 

( + -007\ 

V —.009 / 

Higher-twist 

( +.004\ 

V —.005 / 

Renormalization Scheme 0 

±.001 

Order ai 

±.0003 

Total Theory Error 

±.005 


TABLE III. The total GLS integral and a a for each bin in 
Q 2 . The errors on the GLS are F(stat) ± ( syst ). The errors 
on a s (Q 2 ) are F(stat) ± (syst) ± (thry). Systematic errors 
are correlated in Q 2 . 


(Q 2 ) 

GLS(Q 2 ) 



Q : 

’) 



1.26 

2.39 

± 

.08 

± 

.14 

0.330 

± .023 

± 

042 

± 

.050 

2.00 

2.49 

± 

.08 

± 

10 

0.303 

± .020 

± 

026 

± 

.036 

3.16 

2.55 

± 

.06 

± 

10 

0.287 

± .008 

± 

034 

± 

.026 

5.01 

2.78 

± 

.06 

± 

.19 

0.165 

± .033 

± 

144 

± 

.024 

7.94 

2.82 

± 

.07 

± 

.19 

0.145 

± .061 

± 

136 

± 

.022 

12.59 

2.80 

± 

.13 

± 

.18 

0.164 

± .068 

± 

101 

± 

.014 
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